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ABSTRACT 

The fluctuations in the cosmic microwave background radiation may contain 
deviations from gaussian statistics which would be reflected in a nonzero value 
of three-point correlation function of AT. However, any potential observation of 
the three-point function is limited by cosmic variance, even if a whole-sky map 
of AT is available. Here I derive a general formula for the cosmic variance of the 
three-point function in terms of integrals over the two-point function. This formula 
can be applied to any cosmological model and to any experimental measurement 
strategy. It gives a fundamental lower limit on the magnitude of a measurable 
three-point correlation function, and hence on the measurable amount of skewness 
in the distribution of AT. 
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Now that anisotropies in the cosmic microwave background radiation have been discov- 
ered on large angular scales (Smoot et al. 1992), attention can be focused on the detailed 
properties of these fluctuations. A question of particular importance is whether or not 
their distribution is gaussian. A gaussian distribution has a vanishing three-point corre- 
lation function. An apparently simple thing to check is whether or not this is true of the 
observed anisotropies (to within the experimental error bars). 

However, theoretical predictions of properties of the CMBR are always probabilistic in 
nature, and fundamentally limited by the "cosmic variance" which arises from our inability 
to make measurements in more than one universe (e.g., Abbott & Wise 1984b; Scaramella 
& Vittorio 1990, 1991; Cayon et al. 1991; White et al. 1993). The prediction of a vanishing 
three-point function is no exception, and so comes with a purely theoretical error bar. 

In this letter I address this problem. Assuming that the fluctuations are in fact gaus- 
sian, I compute the expected deviation from zero of the three-point correlation function 
(which includes, as a special case, the skewness), due solely to the effect of making measure- 
ments in a single universe. The answer is expressed in terms of the theoretical two-point 
correlation function, which depends on both the window function of the experiment and 
the cosmological model which is adopted. I give detailed results for the COBE window 
function, assuming a scale invariant (n = 1) power spectrum for the fluctuations. Any 
measured three-point function must have a magnitude exceeding the values given here in 
order to be construed as evidence for non-gaussian fluctuations. I compare my results for 
the cosmic variance of the skewness with previous Monte Carlo simulations (Scaramella & 
Vittorio 1991). 

To begin, let T(n) denote the temperature difference which the experiment assigns to 
a point on the sky specified by the unit direction vector ii, including the effects due to 
finite beam width and any "chopping" strategy which the experiment uses. We can make 
a multipole expansion of the form 

f(n) = a lm W lm Y lm (n) , (1) 
lm 

where the Yj TO (n) are spherical harmonics, the Wi m represent the window function of the 
experiment, and the a; TO are random variables whose distribution must be specified by a 
specific cosmological model. In general, rotation invariance implies that 

( a lm a l'm>) = C l S W S mm> , (2) 

where the angle brackets denote an ensemble average over the probability distribution for 
the a\ m . Assuming an n = 1 power spectrum and ignoring the possible contribution of 
tensor modes, the Sachs-Wolfe effect results in (Peebles 1982; Abbott & Wise 1984a; Bond 
& Efstathiou 1987) 

47 = 1(1 + 1) { > 



2 



for / > 2, where Q is the r.m.s. ensemble average of the quadrupole moment (which may 
differ significantly from the actual, measured quadrupole moment). Since the monopole 
term must be removed from T(n) , we have Wqq = and 

/^f(n) = 0, (4) 

where dQ denotes integration over the unit vector n. The theoretical two-point correlation 
function is given by 

C 2 (ni,n 2 ) = (f(ni)f(n 2 )) 

= J2 C l\ W lm\ 2 Yi m (ni)Y l * m (n 2 ). (5) 
lm 

If Wi m is independent of m, we can simplify this to 

C 2 (ni, n 2 ) = C(x) = i- ^{21 + 1) C x |W,| 2 P^x) , (6) 

where x = n\-h 2 is the cosine of the beam-separation angle and Pi(x) is a Legendre 
polynomial. 

Given the observed values of T(h) for each n, the experimental value of the skewness 
5 is given by the sky average of T 3 : 

s = /§*»(*). m 

(The effect of partial sky coverage will be discussed later.) The corresponding dimensionless 
skewness parameter is J = S/[J <inT 2 (n)/47r] 3 / 2 , but a discussion of this quantity is 
complicated by the need to treat both its numerator and denominator as random variables. 
I therefore concentrate on S as defined in equation (7). 

If the distribution of the a\ m is gaussian, then the ensemble average of the skewness is 
zero: (S 1 ) = 0. Of course, for our particular sky the experimental value of S is unlikely to 
be exactly zero, even if we completely neglect experimental noise (Scaramella & Vittorio 
1991). In order to tell whether or not a particular measured value of S is significant 
evidence of a departure from gaussian statistics, we need to know the cosmic variance 
of S : assuming that the underlying statistics are indeed gaussian. Since the mean of S 
vanishes, we can write the variance of £ as 

^)=j^'^{T\ tn) f\n 2) ). (8) 

Using standard combinatoric properties of gaussian distributions, we have 

<f 3 (n 1 )T 3 (n 2 )> = 9 {f\h 1 )){f\h 2 )){f{h 1 )f{h 2 )) + 6 <T (fn)T(n 2 )) 3 . (9) 
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If W\ m is independent of m, then ensemble averages are always rotationally invariant, and 
(T 2 (ni)) is independent of hi. From now on, for simplicity, we will assume that this is 
the case. Equation (4) then implies that the first term in equation (9) will vanish after 
integration over h\. Thus we have 

<S 2 > = 6/^Ci( ftl , ft2 ). (10) 

Making use of equation (6) , we can simplify this to 

r+1 

(S 2 ) = 3j dxC s (x) . (11) 



This formula assumes full sky coverage. In the case of partial sky coverage, the analysis 
of Scott et al. (1993) can be applied. Specifically, the integrals in equation (10) should range 
only over the solid angle A which is covered, and each factor of 4tt should be replaced by A. 
Since Ci{n\-, %) is always sharply peaked near h\ = fi2, the net change is an enhancement 
of {S 2 } by a factor of 4it/A, provided that A is large enough to encompass the entire peak, 
and provided that equation (4) still holds, at least approximately. For more details in the 
context of a different but similar calculation, see Scott et al. (1993). 

As a specific example, the COBE group (Smoot et al. 1992) has reported results for the 
two-point correlation function with the monopole, dipole, and quadrupole terms removed, 
and the higher moments weighted with a 7°.5 FWHM beam, which results in W\ m = W\ = 
exp[-^(Z + \) 2 cr 2 ] with a = 3°.2 for / > 2, and W\ = for / < 2. Assuming equation (3) 
for the moments of the temperature distribution gives C(l) = 3.65 Q 2 , and performing the 
integral in equation (11) numerically yields {S 2 } = 1.10 for full sky coverage. Including 
the / = 2 term gives instead C(l) = 4.63 Q 2 and (S 2 ) = 3.44 Q 6 . In either galactic 
latitude cut of \b\ > 20° enhances (S 2 ) by a factor of An/ A = 1/(1 - cos 70°) = 1.52. 

Thus, the skewness of the temperature distribution, as determined from the COBE 
data with \b\ > 20° and the quadrupole removed, would have to be significantly larger 
than (5" 2 ) 1 / 2 = 1.3 Q 3 in order to constitute evidence for non-gaussian fluctuations. Here, 
again, Q is the r.m.s. ensemble average of the quadrupole moment, which must be de- 
termined from the normalization of the full two-point correlation function. Q is a con- 
venient measure of the overall magnitude of the fluctuations; for example, in inflationary 
models, the amplitude of the perturbations at the time of horizon crossing is given by 
e H = (12/57r) 1 / 2 <5/To (Abbott & Wise 1984a,b). If we assume an n = 1 spectrum, then 
the COBE data indicates that Q = 16.7 ± 4/jK (Smoot et al. 1992). 

Although the the probability distribution for S must be symmetric about S = 0, it will 
not be gaussian. (It is easy to check, for example, that (S 14 ) does not equal 3(S' 2 ) 2 , as it 
would for a gaussian distribution.) Thus, ±(S' 2 ) 1 / 2 does not correspond to a true 68% con- 
fidence interval. The probability distribution for the dimensionless skewness parameter J 
was computed via Monte Carlo methods by Scaramella & Vittorio (1991). For an n = 1 
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spectrum, beam width a = 3°.0, and including the / = 2 term, they found the J distribu- 
tion to be well approximated by a gaussian with (J7" 2 ) 1 / 2 = 0.2. For comparison, we can 
compute (S 2 ) l / 2 /C{lf/ 2 . This is a ratio of averages rather than the average of a ratio, 
so exact agreement is not expected. Nevertheless, in this case I find C(l) = 4.78 Q 2 and 
(S 2 ) = 3.56 Q 6 , so that (S 2 ) 1 / 2 /C(l) 3 / 2 = 0.18, in good agreement with Scaramella & 
Vittorio's value of (J 2 ) 1 / 2 . 

Experiments on small angular scales typically result in a theoretical two-point cor- 
relation function which can be adequately approximated by C(x) = Coexp[(x — l)/6> 2 ], 
where Cq and 9 C are computable functions of the cosmological model and the experimental 
parameters (e.g., Bond & Efstathiou 1987). In this case, combining equation (11) with the 
correction factor of Scott et al. (1993) for partial sky coverage yields 



[S<) = (4n/A)9 2 C$ 



(12) 



Equation (12) applies when A ^> 9 2 ; the upper limit on (S 2 ) is 6CfJ. 

These results for the cosmic variance of the skewness can be extended to the full 
three-point correlation function. Experimentally, the three-point correlation function is 
determined by first choosing a fixed configuration of three direction vectors; the simplest 
choice is an equilateral triangle, 



77-1 -77-2 = 77.2-773 = 773-771 = cosct , 



(13) 



and we specialize to this case from here on. The experimental three-point correlation 
function is then 

((a) = j dR f(Rh 1 )f(Rh 2 )f(Rh 3 ) . (14) 

Here the hi are to be held fixed in a configuration obeying equation (13), and R is a rotation 
matrix; appropriately integrating over R results in an average over all possible triangles 
obeying equation (13). Note that £(0) is simply the skewness S. A specific realization of 



the hj is 



al 1 



hi = (s a , 



Cry) 

2 
V3 



where Cq, = [(1+2 cos a)/3] 1//2 and s a = [(2 — 2 cosct)/3] This is an equilateral triangle 
centered on the z axis. A specific realization of R is 

f cos <fi sin \ / 1 \ f cos ip sin ip 0\ 



c a) 1 
1/2 



(15) 



R 



sine 
— sin (p cos « 






1/ 



cos 9 sin 9 
V0 -sin# cos 9 J 



— sin ip cos ip 
1/ 



(16) 



This is a product of a sign (e = ±1) and three rotations, the first about the z axis by ip, the 
second about the x axis by 9, and the third about the z axis by <j>. The two possible values 
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of e must be summed over in order to get triangles of both orientations. The integration 
measure in equation (14) for this realization of R is 




Assuming a gaussian distribution, the ensemble average of the three-point correlation 
function vanishes: (C(«)) = 0. The cosmic variance of ((a) is then 

(C 2 (a)) = j dR dR' (f(Rhi)f(Rh 2 )f(Rh 3 )f(R'hi)f(R'h 2 )f(R'h 3 )) . (18) 

Using the rotation invariance of the ensemble average, the integral over R' can be removed, 
leaving 

(CV)} = J dR (T(i?ni)f (Rh 2 )T(Rn 3 )T (ni)f (n 2 )f (%)> . (19) 

Using the appropriate generalization of equation (9), symmetries of the equilateral triangle, 
and / dR T(Rh) = [which follows from equation (4)], we ultimately get 

(CV)) = 6 J dR C 2 (Rn 1 ,n 1 )C 2 (Rn2,n2)C 2 (Rn 3 ,n 3 ) . (20) 

When a = 0, equation (20) reduces to equation (11). 

The value of ((^(ct)) 1 / 2 can be computed by doing the integral in equation (20) numer- 
ically. The result is plotted in Figure 1, assuming an n — 1 power spectrum, the COBE 
window function (7°5 FWHM beam, / = 0,1,2 terms removed), and full sky coverage. 
This indicates the magnitude of ((a) which is to be expected in the COBE data due to 
cosmic variance, even though the ensemble average of ((a) is (assumed to be) zero. The 
experimental value of ((a) would have to rise significantly above the curve in Figure 1 in 
order to indicate a departure from gaussian statistics. The restriction to |6| > 20° in the 
COBE data would enhance ((^(ct)) 1 / 2 by a factor of 1.23 for small values of a, and by 
a somewhat larger factor (which is difficult to estimate reliably) for a > 20°. Note again 
that while the probability distribution for ((a) must be symmetric about ((a) = 0, it will 
not be gaussian, and so ±(C 2 (a)) 1 / 2 does not represent a true 68% confidence interval. 

As an example of the utility of this result, let us examine the effect of cosmic variance 
on the prediction of inflationary models for the three-point correlation function. For an 
n = 1 power spectrum, the angular dependence of the three-point function is determined, 
but the overall amplitude is model dependent (Falk et al. 1993). In single-field inflation 
models, the amplitude of the three-point function is proportional to /u/if, where \x is the 
coefficient of the cubic term in the potential of the inflaton field, and H is the Hubble 
constant during inflation. In two-field models (e.g., Kofman et al. 1991), this amplitude 
is proportional to the strength of the coupling between the two fields. Figure 2 shows 
the ensemble average of the three-point correlation function, assuming the COBE window 
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function (Falk et al. 1993). The overall amplitude is that of a single-field model with 
IxjH = 0.1. The shaded area indicates a band of i^^ct)) 1 / 2 about the mean, assuming 
full sky coverage. Of course, there will be corrections to (£ 2 (a)) coming from the non- 

9/9 

gaussian aspect of the distribution, but these will be suppressed by a factor of /j, /H , 
and cannot significantly change the result. We see that unless fx/H is larger by at least a 
factor of two, the wiggles at large a drop below the uncertainty due to cosmic variance. 
This is unfortunate; Luo & Schramm (1993) have pointed out that the height of the a = 
peak depends sensitively on the power-spectrum index n, but to determine n from the 
three-point function requires normalizing it at large a. This is impossible unless the three- 
point function is large enough to raise its a > 20° features well above the uncertainty due 
to cosmic variance (to say nothing of the uncertainty due to experimental noise, which is 
likely to be at least comparable). Furthermore we see that we must have fi/H > 0.015 in 
order to keep even the a = peak out of the cosmic noise. 

To summarize, the fluctuations in the cosmic microwave background radiation may 
contain deviations from gaussian statistics which would be reflected in a nonzero value of 
three-point correlation function of AT. Any potential observation of the three-point cor- 
relation function is limited by cosmic variance, even if a whole-sky map of AT is available. 
The cosmic variance of the three-point correlation function depends on both the under- 
lying cosmo logical model and the experimental measurement strategy. For COBE, the 
skewness, defined as the sky average of (AT) 3 with the dipole and quadrupole terms re- 
moved from AT, must exceed 1.3 Q 3 in order to be indicative of non-gaussian fluctuations. 

For small-scale experiments, the skewness must exceed (47t/A) 1 / 2 6' c Cq^ 2 . 

I thank Douglas Scott and Martin White for helpful comments. This work was sup- 
ported in part by NSF Grant Nos. PHY-91-16964 and AST-91-20005. 
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FIGURE CAPTIONS 

Fig. 1. — The square root of the cosmic variance of the the three-point correlation function 
((^(a)) 1 / 2 as a function of the beam-separation angle a, in units of (where Q is the 
r.m.s. ensemble average of the quadrupole moment), assuming an n = 1 spectrum of 
fluctuations, the COBE window function (7°5 FWHM beam, I = 0, 1, 2 terms removed), 
and full sky coverage. 

Fig. 2. — The prediction of inflation for the three-point correlation function (((a)) as a 
function of the beam-separation angle a, in units of Q , assuming an n = 1 spectrum 
of fluctuations and the COBE window function. The overall amplitude is that of a 
single-field inflation model specified by fi/H = 0.1. The gray band indicates a range of 
±(C 2 (a)) 1 / 2 about the central value, assuming full sky coverage. 
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